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Size-dependence of Strong-Coupling Between Nanomagnets and Photonic Cavities 

6. O. Soykal and M. E. Flatte 
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University of Iowa, Iowa City, IA 52242 

The coherent dynamics of a coupled photonic cavity and a nanomagnet is explored as a function 
of nanomagnet size. For sufficiently strong coupling eigenstates involving highly entangled photon 
and spin states are found, which can be combined to create coherent states. As the size of the 
nanomagnet increases its coupling to the photonic mode also monotonically increases, as well as the 
number of photon and spin states involved in the system's eigenstates. For small nanomagnets the 
crystalline anisotropy of the magnet strongly localized the eigenstates in photon and spin number, 
quenching the potential for coherent states. For a sufficiently large nanomagnet the macrospin 
approximation breaks down and different domains of the nanomagnet may couple separately to the 
photonic mode. Thus the optimal nanomagnet size is just below the threshold for failure of the 
macrospin approximation. 
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I. INTRODUCTION 

Coupling between an electromagnetic field and an elec- 
tronic transition in matter, with coupling stronger than 
environmental dissipation, has permitted delicate elec- 
tromagnetic control of electronic states. This control al- 
lows sensitive measurement of unknown environments, 
such as the extensive use of nuclear magnetic resonance^ 
(NMR) as a diagnostic probe, as well as the manipulation 
of quantum information, such as in the demonstration 
of atomic teleportatiorP. Successful efforts in this area 
have tended to progress from the systems most weakly 
coupled to the environment (such as nuclei in NMR) 
to more dissipative systems (such as electron spin res- 
onance in solids, first in insulators, later in metals and 
semiconductors). However, even though the dissipation 
is stronger in solids, the coupling is also stronger, sug- 
gesting the potential for very rapid exchange of quantum 
states between light and matter. Recently several ex- 
amples of strong coupling between a single exciton and 
a single photon in a semiconductor have been demon- 
strated, through the mixing of the exciton and photon in 
photoluminescence^El, through Rabi oscillations between 
exciton and photorP, and through optically-induced spin 
rotatiorPH ( sp i n AC Stark effect). 

Multiply-excited atomic systems coherently interact- 
ing with a photon mode exhibit additional unusual phe- 
nomena, such as superradiance^. Multiply-excited ex- 
citonic systems in solids suffer from decoherence due to 
homogeneous and inhomogeneous linewidths and long- 
range Forster coupling between different excitonic tran- 
sitions. Mitigation of both homogeneous and inhomoge- 
neous linewidths is possible by coupling the excitons to 
each other through the coulomb interaction (excitonic 
condensate 11 ), or indirectly through t he ca vity mode, 
in order to form a polariton condensat e! 12 * 13 !. Excitonic 
condensates associated with finite-energy excitons, how- 
ever, are challenging to generate, and are not found at 
room temperature. Strong coupling in multiply-excited 
systems would therefore benefit from a robust, room- 



temperature, coherent electronic state whose coupling to 
the photonic mode can be made larger than its decoher- 
ence rate. 

Ferromagnets are robust room-temperature many- 
body states that couple directly to light, although the 
magnetic dipole transitions associated with individual 
spins couple more weakly to photons than electric dipole 
transitions (by a factor of the fine structure constant^]) . 
Recently it has been pointed oulP^ that the coherent ex- 
citation of the ground-state spin of a small ferromagnet 
(a nanomagnet) can be described by a coupling strength 
orders of magnitude stronger than that of a single ex- 
citonic transition. The locking of the large number 
of constituent spins by the exchange interaction into a 
macrospin causes an increase in coupling strength pro- 
portional to the square root of the number of exchange- 
locked spins. 

Here we expand on the description in Ref. 1151 treat- 
ing with particular care the dependence on nanomagnet 
size of the coupling strength between the nanomagnet 
and the photonic cavity. We find that, for a specific 
magnetic material, the coupling strength increases ac- 
cording to the square root of the volume of the nano- 
magnet (corresponding to the square root of the total 
nanomagnet spin) in the absence of any photons inside 
the cavity. However, when the system is driven in the 
superradiance regime, this coupling strength becomes 
proportional to the volume (or total spin) to the 3/2 
power. We provide estimates of the coupling strength 
for nanomagnets in a spherical cavity and compare with 
the effect of crystalline magnetic anisotropy. The cou- 
pling strengths found are large enough to establish eigen- 
states involving large numbers of entangled photons and 
spin orientation states. These states can be combined 
to generate coherent oscillations of the spin orientation 
and photon number. However, for small nanomagnets 
the crystalline magnetic anisotropy greatly exceeds the 
nanomagnet-cavity coupling, quenching these coherent 
oscillations. For large nanomagnets the macrospin ap- 
proximation, assumed here, fails and the multiple do- 



mains of the nanomagnet separately couple to the cavity. 
The effect of using plasmonic techniques to enhance the 
magnetic field associated with the photonic mode near 
the nanomagnet is described, which may lead to sub- 
microsecond oscillation times for coherent multiphoton 
oscillations in the cavity. 

We begin by describing the nanomagnet-cavity system, 
quantizing the photons of the spherical cavity, and de- 
riving the Hamiltonian of the system. We solve for the 
eigenstates of the coupled system by mapping the dis- 
crete system onto a continuum representation similar to 
a one-dimensional tight-binding model with a spatially- 
varying effective mass. Perturbations to the magnetic 
system such as magnetic anisotropy can be described 
as spatially-varying potentials for this one-dimensional 
tight-binding model. The time-evolution of coherent 
states is evaluated, and the source of dephasing discussed. 



II. 



COUPLED NANOMAGNET-CAVITY 
FORMALISM 

A. Nanomagnet Properties 




FIG. 1: (color online) Schematic of the nanomagnet-cavity 
system with a spherical nanomagnet of radius ro placed at a 
distance of d from the center of a spherical photonic cavity of 
radius R. The orientations of the electric E and magnetic field 
H at the nanomagnet site are shown for transverse magnetic 
(TM) and transverse electric (TE) modes of the photonic cav- 
ity. A uniform magnetic field, Bo, applied along the z-axis 
causes precession of the nanomagnet macrospin, S, with fre- 
quency cj, in resonance with the TM mode of the cavity. 

As shown schematically in Fig. [TJ the oscillator is a 
spherical nanomagnet with a radius r$ possessing a very 
large exchange-locked spin S placed a distance d from 
the center of the cavity for more efficient coupling to the 



cavity mode. Precession of the nanomagnet macrospin 
at a frequency w resonant with the cavity is achieved by 
applying a uniform magnetic field B along the z-axis of 
the cavity. 

A nanomagnet acting as a macrospin, as seen experi- 
mentally in nanomagnet oscillators of roughly this size^, 
has a magnetization 



M = n/V 
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in terms of the collective spin operator S and the Heavy- 
side step function Q(x). The magnetization in Eq. (fTl) 
depends on the spin density of the nanomagnet. The 
modal coupling (the coupling of the nanomagnet to the 
photonic mode) is the overlap of this magnetization with 
the cavity mode amplitude. For a nanomagnet that is 
small in size compared to the length scale of variations 
in the cavity mode strength, the coupling will be inde- 
pendent of the spin density and will only depend on the 
total spin. It is possible, however, to enhance the modal 
coupling through mode design, such as is common to en- 
hance the interaction between gain media and an optical 
cavity in semiconductor lasers^. For example, an optical 
field is stronly enhanced near a sharp metal object (used 
in tip-enhanced spectroscopy^); a similar approach here 
could be used to strongly enhance the strength of the 
nanomagnet-cavity coupling. 



B. Quantized electromagnetic field in a spherical 
cavity with a nanomagnet 

The presence of the nanomagnet in the cavity, and its 
magnetization field, modifies the properties of the dy- 
namic electromagnetic field in the cavity. The nanomag- 
net precesses in the static external magnetic field, yield- 
ing a temporally-oscillating magnetization characterized 
by the precession frequency uj. Thus the nanomagnet be- 
haves as an oscillating source in the Maxwell equations 



V H = 0, 

V E = 0, 



V x E - iky/^ /e Q H = 0, 
V x H + ikE/^no/eo = V x M. 

(2) 
We introduce the time dependence of the fields (e 1 "*) into 

H, E, and M . This produces the following Helmholtz 

wave equations, 



k 2 ) (r ■ H) 



-%L ■ (V x M) 



k 2 ) (r-E) = Z kL ■ M. 



(3) 



From these the solutions of the transverse magnetic (TM) 
and electric modes (TE) can be obtained 
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where the vector spherical harmonics Yi^^ n are defined 
as LYi m (8, (f>)/ \Jl(l + 1), in terms of angular momentum 
operator of the field L and spherical harmonics Yi m \^ln 
the most general form, they are defined as 



Yj,i, mj = C(l, 1; mi, m\j, mj)Yi r , 



(5) 



in terms of Glebsch-Gordan coefficients and helicity basis 
vectors e m . The helicity basis vectors form a spherical 
tensor of rank 1, i.e. 



(6) 



e ± = T(x±iy)/V2, 



where eo = z. 

The fi(kr) and gi(kr) appearing in Eq. Q are the 
corresponding solutions for the radial part of each mode, 
A\ ' (kr)h\ (kr)+A\ ' (kr)h\ (kr), in terms of the spher- 
ical Hankel functions. The coefficients a lm and a\ m , 
which specify the amounts of transverse magnetic and 
transverse electric multipole (l,m) field strengths, are 



ATM) 
x lm 



(TE) 



ik 3 



-k 2 



j l (kr')YZ n (0',<l>')L.Md l r', (7) 



vw 



=? j 3iW)YC m W ,#)L • (V x M) dV, 



which holds for any well-behaved vector field. Therefore, 
the TM mode will be the only non- vanishing mode to be 
considered in our interaction Hamiltonian. The condition 
H± = is trivially satisfied by r ■ U; m = at the cavity 
walls, whereas the condition Eu = gives 

[r x (V x fi(kr)Y ltltm )] \ r=R = -d r (rf l (kr))Y ltl>m = 0. 

(9) 
For waves that are finite at the origin, the suitable 
choice of fi(kr) is the spherical Bessel function of first 
kind, ji(kr). Hence, the normalization integral of the 
basis functions tt/ m = ji(kr)Yij tm , modified for xi y (cor- 
responding to the zeros of \rji(kr)\'), 

/ U*m u l'm'd 3 r = / ji(kr)ji,(kr)Yl lm Yi^i^ m ,r 2 drdn 
= y f 1 - ^p ) \Ji(xh)\ 2 8u'8mm', (10) 



yields to the following mapping of the multipole strength 
coefficients onto cavity photon creation and annihilation 
operators: 



where the volume integration is carried over the local 
sources. 

The radiation of the cavity field is due to the harmoni- 
cally oscillating components of the nanomagnet magneti- 
zation in the x-y plane orthogonal to the radial direction, 
namely M x>y . Because of this specific symmetry of the 
cavity-nanomagnet system, the multipole field strength 
coefficients a, for TE mode will simply vanish due to 
the relation 

L-(VxM) = iV 2 (r-M)--^(r 2 V-M), (8) 
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which satisfy the appropriate Wcyl-Hciscnbcrg commu- 
tation relations, [ai m ,a\, m ,] = 8w5 mm >- 

Therefore, the second quantized form of the magnetic 
field for the cavity TM mode becomes 



h (tm) = \p l 

f^\Ji(x h )\ 



1(1 + 1) 



*h 



-1/2 



hu. 
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VoR 3 



[ a hn U lm + a lm U ^ 



(11) 



The total Hamiltonian of the system incorporates the magnetization M of the nanomagnelP^, 
magnetic H and electric E fields of the cavity and the 

JT = - I (no\H\ 2 + e \E\ 2 +^(H- M)) d 3 r. (12) 



The first two integrands on the right hand side of Eq. ( 12 ) 



correspond to the free field Hamiltonian, whereas the 
third integrand is the interaction Hamiltonian of the 
nanomagnet-cavity system, 



^ = £r 



(TM) (TM) 



'a 



tin 



M ■ ui m d 3 r + c.c. (13) 



with the coupling constant, 



n(TAf) 



2\M*h)\ 



1(1 + 1) 



-1/2 



R 3 



, (14) 



for a TM mode with angular momentum I. All compo- 
nents of the field are identically zero if / = m = 0, a re- 
sult associated with the absence of radiating monopoles. 
From Eq. Q that dipolc field strength coefficient (1 = 1) 

dominates over other multipoles, i.e. a\ m » a 2m 3> 

a 3m -^ ■ ■ ■ ■ The basis functions for the dominant 
dipole TM mode (I = 1) are 

un = ^=h(kr)(Y 11 (0,cf > )e o -Y lo (e, <!>)$+), 



Uio 



un = 
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V2 
-1 
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n (kr) (Y xl (6, 4>)e_ - Y xl (6, 4>)e+) , (15) 
j 1 (kr)(Yu(e,<f>)e -Y 10 (e, <!>)$-), 



C. Coupling of the Nanomagnet to the Photonic 
Cavity 

To describe the coupling of the nanomagnet to the cav- 
ity, the spin operators of the nanomagnet should be writ- 
ten in the same helicity basis as the photonic field, 

S = -^=(S+e_ - S-e + ) + S z e Q , (16) 

in terms of the nanomagnet spin raising and lowering 
operators 



S±\l Sl m s ) = y/(l s =F m s )(l s ± m s + l)\l s ,m s ± 1). (17) 

Introduction of this total spin operator S and the basis 
functions of the spherical wave expansion (Eq. (15)) into 
Eq. (13) yields a fully quantum treatment of the total 



Hamiltonian for the nanomagnet-cavity system, 



J#L ~ Hlu~ [ aLa-, + 



1 



1 ■ • 2J ' " n 

gHB~T~{ (a 7 S + + a^S-) 



g—B s z 



(18) 



in which the spin interacts only with a single photon 
mode 7. Modes of higher I would be out of reso- 
nance because of the cavity quantization, and energy 
non-conserving terms with negative helicity have been 



dropped (relying on the rotating wave approximatiorPU) . 
The nanomagnet-photon coupling constant T 7 becomes 



Ji(kd) 
8ftb"i(l/i 7 )| 



1(1 + 1) 



-1/2 



3fi.w 7 /io 
ttR 3 ' 



(19) 



where the mode frequency w 7 is related to the radius of 
the cavity R with fci 7 = wi 7 /c = Xij/R. 

The interaction with the uniform magnetic field Bo, in- 
troduced in Eq. ( 18 ), sets the cavity in resonance with the 



energy level splitting of nanomagnet spin states when- 
ever the relation Sw 7 = g^BB® is satisfied. Therefore, 
any spin flip up (down) process of the nanomagnet spins 
results in an absorption (emission) of a cavity photon in 
the case of exact resonance, e.g. an applied uniform mag- 
netic field of So = 7 T, corresponding to a precession of 
the macrospin with a frequency of ~ 200 GHz, will cause 
the nanomagnet spins to be in exact resonance with a 
cavity volume of 1.25 mm 3 . We assume the lowest TM 
mode of the cavity is in resonance with the spin-flip tran- 
sitions of the nanomagnet, so as higher-energy modes will 
not be in resonance the subscript 7 will be omitted from 
Eq. ([18]). 

The eigenstates of the nanomagnet, treated as a 
macrospin, are simultaneous eigenstates of the total spin 
operators S 2 , and S z given by \l s ,m s ), where \m s \ < 
Is — N/2. Part of the macrospin approximation is the 
assumption that l s is fixed, and most likely it will be the 
maximal spin state l s = N/2 due to additional energy re- 
quirement of any other l s ^ N/2 subspace. The Hilbert 
space of N independent spins should include the states of 
a macrospin corresponding to l s = N/2. Therefore, the 
structure of these basis states is similar to those of the 
Dicke modePSl for iV independent atomic spins, wherein 
l s is the cooperation number of the paramagnetic collec- 
tion of spins. However, for a realistic nanomagnet, ele- 
ments of the Hilbert space with l s 7^ N/2 are split off in 
energy due to the exchange interaction giving rise to ex- 
tra mechanisms, i.e. elementary excitation of spin waves 
(magnons) . 

Since each magnon excitation reduces the total mag- 
netic moment (fi o< l s ) of the nanomagnet in the amount 
of 2.21/ib for Fe, it is possible for the nanomagnet to- 
tal spin angular momentum to start in a different l s 
subspace rather than the maximal l s = N/2. This re- 
duction in Z s is less than 1% at room temperature for 
iron, suggesting that nanomagnet oscillators of approx- 
imately these sizes can be well-described by as having 
maximal spin at room temperature. The validity of the 
macrospin approximation relies on the effectiveness of the 
exchange-locking of the spins at room temperature. For 
the nanomagnets we consider here, spherical nanomag- 
nets of radius r$ = 2.3 nm, 11 nm, and 50 nm consisting 
of iron (magnetic moment 2.21/zs per atom), and pos- 
sessing N ~ 10 4 , 10 6 , and 10 8 electron spins, respectively, 
the macrospin approximation is reasonable^ (although 
perhaps questionable at for the largest nanomagnet con- 
sidered) . 



III. PROPERTIES OF THE COUPLED 
NANOMAGNET-CAVITY HAMILTONIAN 

The total excitation number 2£, corresponding to the 
maximum number of photons n in the cavity (when the 
nanomagnet is parallel to the static magnetic field), needs 
to be conserved by the Hamiltonian in Eq. (m). For an 
initial configuration of the macrospin pointing antiparal- 
lel to the static field Bo and no photons in the cavity, 
£ = N/2, the basis states of the spin-photon mode sys- 
tem \n,m s } can be written as \n, £ — n) or |£ — m s , m s ), so 
that the basis states are indexed either solely by photon 
number of the cavity (n), or by eigenvalue of S z (m s ). 

To proceed, we adopt the notation \n, £ — n) and drop 
the redundant reference to the m s , so the total Hamilto- 
nian takes the form of 

.^ = Y^ Eo\n)(n\-T(n) [\n + l)(n\ + \n)(n + 1|] , (20) 

n=0 

in the Fock space, where the constant energy coefficient 
Eq term and the coupling strength r(n) are defined as 



E = M£ + l/2), 
r(ra) = tiTgii B (n + l)y/2£ 



n . (21) 

In matrix form, the same Hamiltonian can be written as 



Jf 



( E Q - T (0) 
-r(0) E -r(l) 
-r(l) E 



V o 







\ 





-r(2£-l) Ej 



(22) 



similar to the Hamiltonian matrix expected for a nearest- 
neighbor tight-binding model with a spatially-dependent 



mass (see Fig. v\. For 2£ 



N = 10 s , the magnet- 
microwave mode coupling, t(ti), changes over a range 
of 33 kHz - 1.3 THz through all possible photon (spin) 
numbers. 6V(n)/6ro = r'(n) acts like a driving force for a 
fictitious particle moving between sites labeled by photon 

number n, so |0) —> — > \n — 1) — » \n) —>■ \n + 1) —¥ 

— > |2£). The solutions n of r'(n)| no = are equilib- 
rium points in cavity photon number, and for this system 
there is one at no = (4£ — l)/3. The coupling can also 
be expressed in terms of the collective spin number m s 
as t(iti s ) = frTgnBiS, — ms + 1)VC + to «j with an equi- 
librium point of too = (1 — £)/3. For a system consisting 
of a very large number of spins (£ 3> 1), the eigenfunc- 
tions of the Hamiltonian in Eq. ( 10 ) are expected to be 
centered about n = 4£/3 as well as m = — £/3. 

For an initial state |n, m s ), if we are only interested 
in transitions which conserve energy and in which a 
photon is emitted, the rate of photon emission R n is 
proportional to X)v* \(^\ a1< S-\n,m s )\ 2 , where |\t) rep- 
resents the possible final states of the system. There- 
fore, R n = A(n + 1) 2 (2£ — n), or equivalently R n — 
A(£ — to s + 1) 2 (£ + to s ). The factor A can be identified as 



4. M a /ll 



n-2 n-1 11 u+ 1 n+2 



i:(n-2) i<n-l) t;(ii) r(n+l) 

FIG. 2: (color online) Lattice-like schematic of the spin-cavity 
Hamiltonian in Eq. ( 22 1 where successive lattice sites repre- 



sent the possible photon states in the cavity. Note that con- 
servation of total excitation number £ can be seen from the 
addition of arrows belonging to the nanomagnet spin states 
along the z-axis S z (purple, long arrows) and the correspond- 
ing cavity photon number n (red, short arrows) for each site. 
Transitions between successive photon states (lattice sites) 
are governed by the magnet-microwave mode coupling t(ti), 
similar to the hopping in a tight-binding model. 



the Einstein A-coefficient by applying R n to a single spin 
pointing upward (£ = m s = 1/2) when the cavity has no 
photons (n = 0). Since R n reaches its maximum value 
of 4^4(iV/3) 3 for the equilibrium point m (or n ) in the 
large spin limit, the equilibrium points no and too are the 
photon number and spin number, respectively where the 
nanomagnet-cavity system exhibits superradiance^. 



A. Solutions in the Continuum Limit 

For N = 10 4 , 10 6 , and 10 8 the solutions of the 
nanomagnet-cavity Hamiltonian corresponds to the diag- 
onalization of large matrices in the form of Eq. ( 22 ) with 
increasing ranks of 10 4 , 10 6 , and 10 8 for nanomagnets of 
radius r$ ~ 2.3nm, llnm, and 50nm, respectively. The 



magnet-photon coupling strengths at the superradiancc 
regime t(toj) are estimated to be roughly 5.3 neV, 5.3 
fJ.eV, and 5.3 meV for these three different nanomagnet 
sizes. 

The eigenfunctions of the nanomagnet-cavity Hamilto- 
nian given in Eq. ( 20 ) can be expanded 



*,= 



E 



$\ri), 



(23) 



in terms of Fock number states and the respective phase 
constants defined by ip™ . Applying the Hamiltonian in 
Eq. (20) onto these states with the aid of Schrodinger 



equation J^tyj = Ej^j, where Ej are the eigenvalues of 
the nanomagnet-cavity system, yields the following re- 
cursion relation 



(Ej - E )^ + r(n - l)^" 1 + r(n)^ 



, r 



0, (24) 



for the phase constants. Since the nanomagnet posseses a 
very large number of spins, the continuum limit consists 
in making the replacement ip 71 — ¥ ipj(ne) for the discrete 



phase constants in Eq. ( 24 1 . Then a continuous lattice- 



like relation can be obtained, 

Ejtpj(ne) + T(ne)i(jj(ne + e) + r(ne — e)ijjj(ne — e) = 0, 

which can also be transformed into the ordinary differ- 
ential equation 



. .dipj(x) dr(x) dtpj(x) 
t[x)- 



dx 2 



dx dx 



(26) 



, . dr(x) 1 d 2 r(x) „ , 



with boundary conditions ipj(0) — ipj(2s z ) = 0, by 



Taylor-expanding the phase constants tpj in Eq. ( 25 1 up 



to C(e 3 ) and defining ne = x. Some of the lowest-lying 
energy eigenvalues Ej and eigenf unctions ipj(x) of this 
differential equation, shown in Fig. HI can be obtained in 
the WKB approximation from 



S(Ej) 




(27) 



where the effective potential is given by V e (x) — r'(x) — 

2 n * r - 1 

r (x)/4t(x) — 2t(x) (see Fig. pj). Shown in Fig. Ware 

eigenstates of the coupled nanomagnet-cavity system for 

three different sizes of nanomagnet. 
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FIG. 3: (color online) The effective potential of the magnet- 
photon system in the WKB approximation is shown with re- 
spect to cavity photon number n centered about the superra- 
diance regime no. 



IV. NANOMAGNET-CAVITY COHERENT 
DYNAMICS 

A. Form of the Coherent State 

A coherent state for the nanomagnet-cavity system 
can be written as a displaced nanomagnet-cavity ground 
state (by a photon number xq from the equilibrium point 
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FIG. 4: (color online) Wave functions of the nanomagnet- 
cavity system shown as a function of photon number, n, cen- 
tered about no, for nanomagnets of radius ro = 2.3,11,50 
nm, consisting of N = 10 4 , N = 10 6 , and N = 10 8 spins 
respectively. First row (a)-(e)-(i) are the ground states with 
a full width half maximum (FWHM) represented in photon 
numbers, second row (b)-(f)-(j) are the first excited states, 
third row (c)-(g)-(k) are the second excited states, and the 
fourth row (d)-(h)-(l) are the 150" 1 excited states. 



n ), or 



ipo(x) 
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-((x-x )-n a ) 2 /2a 2 



(TV2TT 



(28) 



where the standard deviation a can be found by matching 
the values of the full width at half maximum (FWHM) of 
the ground state and the Gaussian function to each other: 
for instance FWHM [ip (x)] = 2V21n2cr ~ 14609 (Fig. Q 
(i)) for N — 10 8 . The eigenfunctions of the nanomagnet- 
cavity system are complete and orthonormal, hence they 
serve as a suitable basis to expand any coherent state 
over, 



Jo 



^x,t) = T i A j e- iE ^%(x). 



3=0 



(29) 



Equating the Gaussian function in Eq. ( 28 ) to the coher- 



ent state of Eq. (29 1 at initial time t — 0, i.e. 
1 



^Ajtpjix) = 



3=0 



aV2n 



,-((x-xo)-n y/2<T- ! 



(30) 



multiplying both sides by i/j'j and using the orthonormal- 
ity condition of the nanomagnet-cavity wavcfunctions re- 
veals the phase constants Aj of the expansion as 



A, 



1 



CTV27T 



2i 



' iP J (x)e- (ix - x ^- n ^ 2 / 2a2 dx. (31) 



For three sizes of the nanomagnet, the coherent states 
shown in Fig. [5ja)-(c), are characterized by large oscilla- 
tions over ranges of (2x =) 1780, 1.76xI0 4 , and 1.76xl0 5 
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FIG. 5: (color online) (a) Amplitude of a coherent state 
for 3 different nanomagnet-photon systems consisting of (a) 
N = 10 4 , (b) N = fO 6 , and (c) N = fO 8 spins are shown as a 
function of photon number n. The large oscillations of these 
coherent states occur about (a) no ~ 6666 with a period of 
T = 1.5ms, (b) n ~ 6.66 x 10 5 with a period of T = 150/xs, 
and (c) no ~ 6.66 x 10 7 with a period of T — 15/is, respec- 
tively. 



photons with periods of T = 1.5 ms, T = 150 /is, and 
T = 15 /is, respectively. Summation over the first 150 
eigenstates (jo = 150) extracted from WKB is sufficient 
enough to obtain convergence in the dynamical proper- 
ties of these nanomagnets. The Zeeman energy of the 
nanomagnet and transverse magnetic field amplitude of 
the cavity at the nanomagnet 's location can also be eval- 
uated from 



(AE Z ) - (<Kx,t)\n z B \<f>(x,t)), 
(Br) = (4>(x,t)\H TM {d)\<f>(x,t)), 



(32) 



respectively, by using the same coherent state represen- 
tation. Large oscillations of these quantities shown in 
Fig. [6] indicates the coherent energy exchange occuring 
back and forth between photons in the cavity and the 
spin states of the nanomagnets. 
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FIG. 6: (color online) Time evolution of the Zeeman energy 
of the nanomagnets (red, solid) consisting of (a) N = 10 4 , (b) 
N = 10 6 , and (c) N — 10 s spins are shown in coherent state 
representation as well as the amplitude of the transverse mag- 
netic mode of the cavity field (blue, dashed) at nanomagnet 
location z = d. 



B. Dephasing of the Coherent State 

The coherent properties of this nanomagnet-photon 
system will also depend on the dephasing of the coherent 
state 4>( x, t), due to inhomogenity of the coupling r(n) in 
Eq. ( |21[ ). The dephasing time of the nanomagnet-cavity 
coherent state can be extracted by a Gaussian fit to the 
peak values of the autocorrelation function between a co- 
herent state at time t and its initial state at t = 0, 



p(t) = Mx,tmx,o))\ 2 , 

2 



(33) 






2 iEjt/h 



whereas each peak (inset of Fig. [7]) is representing the 
revival amount of the coherent state after every success- 
ful period T of oscillation. Exceptionally long dephas- 
ing time of order seconds are shown in Fig. [7J As the 
nanomagnet gets bigger the change in r(n) with n be- 
comes smoother and smoother, leading to longer dephas- 
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FIG. 7: (color online) Dephasing time of the coherent state 
for nanomagnet-photon systems of (a) TV = 10 4 , (b) TV = 10 6 , 
and (c) TV = 10 8 spins (or equivalently photons) obtained by 
a Gaussian fit to the peak values of the dephasing functions 
(insets) at successive time intervals. Each peak value repre- 
sents the amount of correlation after every full period T of 
oscillation. 



ing times. 

Although this treatment is for zero temperature, the 
coherent properties of the nanomagnet-photon system 
should persist to as high a temperature (and over as long 
a timescale) as the macrospin description remains reli- 
able. We have assumed an infinite Q for the cavity, so the 
decoherence of the system is expected to be determined 
by photon leakage from the cavity, rather than these ex- 
ceptionally long calculated times. Furthermore, the el- 
ementary spin excitations (magnons) would not directly 
affect the dephasing of the system, for magnons preserve 
the spin quantum number m s , requiring an up spin to flip 
down for every down spin flipping up. In realistic nano- 
magnets, spin-lattice coupling of m s to phonons through 
spin-orbit coupling will cause a cutoff of the dephasing 



times shown in Fig. [7| For spheres of yttrium iron garnet 
(YIG) at low temperature this spin-lattice time is sev- 
eral /us! 23 * 24 ! Therefore, observation of a full oscillation 
cycle should be possible for nanomagnets with a radius 
of 50 nm or larger. On the other hand, the times at room 
temperature in YIG (~ 200) ns^and iron (~ 20 ns)P^are 
too small to observe a full oscillation. However, coherent 
dynamics corresponding to a portion of the oscillation 
involving ~ 24 photons/ns, or ~ 470 photons for iron 
and 4700 photons for YIG should be still observable for 
the nanomagnet with radius r = 50 nm. If, however, the 
modal coupling is increased using approaches such as tip- 
enhancement of the optical field, then the coupling could 
be far stronger even for a small nanomagnet. Guided by 
estimates from tip-enhanced Raman spectroscopy 18 , the 
intensity of the mode at the nanomagnet 's position could 
be increased by 10 2 — 10 6 , leading to enhancements of the 
oscillation frequency of order 10 — 10 3 . 



C. Crystalline Magnetic Anisotropy 

We also examine other deviations from ideality for the 
nanomagnet, such as the spin dependent cubic crystalline 
magnetic anisotropy (CMA). The CMA of iron is given 
by 



E, 



CMA 



U\ (k 2 K2 + K 2 K 3 + K \ K V) + U2K1K2K3, (34) 



where U\ — 4.2 x 10 5 erg/cm 3 and XJ% — 1.5 x 10 5 erg/cm 3 
are the cubic anisotropy constants for iron at room tem- 
perature and an arbitrary magnetization direction is de- 
fined by the directional cosines «i,«a,«3 referred to the 
cube edges. Since the nanomagnet is a sphere, shape 
anisotropy is not relevant. In the case of a cubic crys- 
tal whose easy axis is aligned along the body diagonal, 
Ecma energy depends on the orientation of the nano- 
magnet spin S, defined by Kj. 

The CMA of iron causes a detuning of the energy spac- 
ing for different spin orientations from the resonant fre- 
quency of the cavity, along with a dispersion in that spac- 
ing. The uniform detuning, corresponding to a uniform 
shift in the precession frequency of the nanomagnet, can 
be compensated for with a slight adjustment in the ap- 
plied magnetic field. The dispersion, however causes a 
variable detuning of roughly 200 neV, 13 neV, and 1.3 
neV of the E$ in Eq. (20) over the range of oscillation 
shown in Fig. [5ja)-(c), respectively. For the smallest 
nanomagnets the effect of CMA dominates over the cou- 
pling between the photons and the spin. For example, 
for a nanomagnet radius of 2 nm consisting of 10 4 total 
spins, the CMA is significantly larger than the magnet- 
photon coupling strength r(n) (~ 5.3 neV) in Eq. (20). 



Therefore the CMA will cause the eigenstates to localize 
in photon and spin number, producing rapid decoher- 
ence for a coherent state. We note that this observation 
largely rules out the possibility of observing these coher- 
ent oscillations in a single molecular magnet^Sl, for the 
spins of these molecules are considerably smaller than 



the spin of the nanomagnet considered above. However, 
this detuning is much smaller than the magnet-photon 
coupling strength of other nanomagnet sizes (10 nm and 
50 nm in radii) and therefore will not destroy the coher- 
ent oscillations for them, although it may still limit the 
dephasing times to shorter than that shown in Fig. ff\h)- 
(c). 



V. CONCLUDING REMARKS 

Calculations for three different nanomagnet sizes in a 
photonic cavity indicate that strong-field coupling be- 
tween photons and spins is possible, and should sub- 
stantially exceed the coupling observed in solids between 
orbital transitions and light. The Hamiltonian for the 
coupled nanomagnet-cavity system is solved in the con- 
tinuum limit to obtain a coherent state representation 
of the system around the superradiance regime. This 
coherent state is characterized by large oscillations in 
photon number of the cavity (or equivalcntly the total 
spin number of the nanomagnet) with exceptionally long 
dephasing times and is expected to be observable for re- 
alistic nanomagnets with radii from 10 — 50 nm. Ap- 
proaches to enhance the coupling, such as using a metal 
tip to enhance the optical field, have been proposed. For 
the smallest nanomagnet (2 nm radius) the dispersion 
caused by crystalline magnetic anisotropy would largely 



quench the coherent oscillations, but for nanomagnets in 
the 10 — 50 nm radius range the coupling to the cavity 
is much stronger than the dispersion caused by CMA. 
The dephasing times increase with increasing nanomag- 
net size, due to the greater uniformity of the coupling 
terms between states that differ by one photon and one 
spin flip. Thus the most coherent nanomagnet-cavity sys- 
tems will be those that are just under the size threshold 
where the macrospin approximation ceases to be accu- 
rate. The effects of magnons have been considered and 
shown to not substantially modify these results. 

Future work shall investigate how to use the strong 
coupling features described here to transfer coherently 
states of the electronic system to the photonic one, and 
back again. A particularly interesting direction will be 
to consider the effect of active nanomagnetic systems, 
such as those demons trate d to be coherently driven by 
electrical spin current j- 16 l 27 H 3]] on ^ ne optical state of the 
cavity. As phase-locking has been demonstrated between 
two such oscillator s] 32 1 33 1 , mediated perhaps by spin waves, 
would it be possible to phase-lock them through interac- 
tion with a cavity such as the one considered here? 
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